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APPLICATION  OF  THE  THEORY  OF  SIMILARITY  TO 
TO  THE  EXPERIMENTAL  INVESTIGATION  OF  THE 
WORK  REGIMES  OF  GAS  TURBINES 

G.  E.  Kalinin 

Similarity  conditions  of  gas  turbines  are  examined; 
similarity  criteria  for  their  work  regimes  on  dif¬ 
ferent  working  substances  are  proposed. 

At  present  the  theory  of  similarity  is  rightly  winning  a  larger 
and  larger  place  in  the  investigation  of  turbine  engines  as  the  sci¬ 
entific  basis  of  the  organization  of  the  experiment.  Based  on  the 
theory  of  similarity,  criterion  relations  have  been  worked  out  for 
generalizing  the  results  of  tests  and  extending  these  results  to 
turbine  work-regimes  not  experimentally  verified. 

In  some  cases,  however,  when  determining  the  indices  of  turbine 
performance  in  accordance  with  the  results  of  testing  it  on  a  working 
substance  different  from  the  natural  one,  the  theory  of  similarity 
cannot  be  used  for  investigating  the  turbines  because  of  the  lack  of 
criterion  equations  for  simulation  on  various  working  substances.  This 
last  is  connected  with  the  fact  that  of  five  criteria  of  similarity 
at  present  established  and  determining  the  turbine  work-regime  only 
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two  may  be  used  for  the  graphic  representation  of  characteristic  tur¬ 
bines.  The  remaining  three,  namely,  the  criteria  of  Reynolds  (Re) 
and  Prandtl  (Pr)  and  the  index  of  the  adiabatic  curve  (k)  must  be 
the  same  in  nature  and  in  the  model ;  but  if  any  one  of  them  differs, 
the  independence  in  practice  of  the  examined  process  from  that  criterion 
must  be  demonstrated.  Values  determining  the  physical  properties  of 
the  working  substance  (viscosity,  heat  conductivity,  compressibility) 
enter  the  criteria  Re,  Pr,  and  k.  Consequently  the  demand  for  equality 
between  the  criteria  Re,  Pr,  and  k  in  nature  and  in  the  model  means 
that  the  turbines,  to  which  are  applied  the  methods,  based  on  the 
theory  of  similarity,  of  generalizing  and  extending  the  experimental 
data,  must  work  on  the  same  working  substances.  In  addition,  that 
means  that  the  temperature  and  pressure  of  the  working  substance  must 
be  close  in  nature  and  in  the  model,  since  the  physical  properties  of 
a  gas  change  with  a  change  in  thermal  parameters.  For  simulation  on 
various  working  substances  it  must  be  possible  to  take  Re  ^  idem,  pr  ^ 
idem,  k ^  idem,  i.e.,  assume  the  independence  of  the  process  from 
these  criteria. 

If  in  fact  the  process  is  not  related  to  any  criterion  it  is 
called  auto-simulating  with  respect  to  this  criterion.  At  present 
it  has  been  ascertained  that  when  Pr  2  5  •  10s  to  7  •  105  the  process 
in  the  turbine  may  be  considered  auto- Simula ting  with  respect  to  the 
Reynolds  number  [1],  There  are  methods  of  accounting  for  the  Influence 
of  the  Re  numbers  on  the  characteristics  of  the  turbine  when  Re  <  5  •  IQ5. 
In  cases  met  in  practice  the  flow  of  gas  in  the  flow-through  part  of 
turbine,  expressed  in  the  criteria  "used  at  present,  cannot  be  con¬ 
sidered  auto- simulating  also  with  respect  to  k.  A  change  in  the  mag¬ 
nitude  of  k  exerts  an  influence  on  the  characteristics  of  the  turbine 


FTD-TT-62-I456/I+2+4 


-2- 


in  all  regimes.  The  matter  of  simulation  when  k  ^idem,  notwithstanding 
research  done  in  this  area,  is  not  resolved  and  is  an  obstacle  to 
simulation  on  various  working  substance. 

In  approaching  the  solution  of  this  problem  it  is  advisable  at 
first  to  examine  the  conditions  of  similarity  of  work  regimes  of  gas- 
turbine  engines.  The  following  may  be  referred  to  the  equation  defin¬ 
ing  the  work  of  the  gas  in  the  turbine:  gas-dynamic  equations  for 
compressible,  viscous,  and  heat-conductive  liquids  [2];  energy  and 
discontinuity  equations  connecting  the  parameters  of  the  gas  on  enter¬ 
ing  and  leaving  the  turbine  with  the  useful  work  taken  off  the  turbine 
rotor. 

The  gas-dynamics  equations: 


X  o  L 

8  '  d'  I 


2  <K 


i _ f _X  • 

*dx9\jdxaJ  g  ’dx?’ 


pv~RT; 


the  energy  and  discontinuity  equations: 


_i 
A 

Yo^o^oa  =  Ya^i^j 


2  a> 


n 


where c  is  the  absolute  velocity  of  the  gas  (the  index  a  designates  the 
component  of  velocity  parallel  to  the  axis  of  the  turbine);  7)  specific 
weight;  g)  acceleration  of  the  force  of  gravity;  t)  time;  jj.)  first 
coefficient  of  viscosity,  taking  displacement  deformations  into  account; 
£)  second  coefficient  of  viscosity,  taking  expansion  deformations  into 
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account;  f)  potential  of  the  force  field;  p)  pressure;  v)  specific 
volume;  R)  gas  constant;  T)  absolute  temperature;  L  )  work  at  the  rotor 
axle;  A)  heat-equivalent  of  the  physical  work;  i)  heat  content;  P) 
area  of  a  cross  section;  x^,  x^,  x^)  coordinates  of  a  point  in  space 
(three-dimensional  space  is  under  consideration,  therefore  p,  v,  and 
a  pass  through  the  values  l,  2,  and  j5;  corresponding  to  this  the 
subscripts  £,  v,  a  with  vectors  will  assume  the  values  1,  2,  and  j5)  . 

Instead  of  the  energy  equation  the  following  equation  may  also 
be  used: 

A  •  75AT , = %iG  (/'*  ~i2'), 

where  N^)  internal  power  of  the  turbine;  G)  consumption  of  gas;  nol) 
internal  efficiency  of  the  turbine;  iD)  full  heat  content  of  the  gas 
before  the  turbine;  ±|)  heat  content  at  the  end  of  isentropic  expansion. 

The  following  may  be  referred  to  single-valued  conditions  of  the 
process  in  the  turbine  (a  turbine  work-regime  which  has  established 
itself  is  under  consideration,  I.e.,  a  process  without  initial  con¬ 
ditions)  :  physical  magnitudes  characterizing  the  working  substance, 
pressure  and  temperature  of  the  working  substance  entering  the  turbine 
and  pressure  on  leaving  the  turbine,  the  geometrical  dimensions  of  the 
flow-through  part  of  the  turbine,  the  angles  of  entry  of  the  stream 
into  the  vanes  of  the  runner  of  the  turbine  (the  angles  between  the 
vector  of  relative  stream  velocity  in  front  of  the  runner  and  the 
vector  of  the  circumferential  velocity  of  the  working  wheel) . 

As  a  result  of  investigating  the  invariance  to  similar  trans¬ 
formations  of  the  enumerated  equations,  of  examining  the  conditions 
of  uniqueness,  of  using  the  method  of  substitution  in  forming  criteria, 
and  of  excluding  constant  criteria,  we  may  note  the  following  criteria 
determining  the  work  regime  (when  auto- simulation  with  respect  to  Re 
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and  Pr  Is  assumed) 


U 1  U ‘  p,  k 

kgRT^kRT,'  p,  ’  T^\  ' 


(1) 


where  u  is  the  circumferential  velocity;  and  the  ascertainable  criteria 
are 

A75N,/  Nt\ 

Gi,  \  0r  QiJ'  1 ol  etS. 

lr 

In  the  quoted  system  the  criterion  does  not  differ  in  prin¬ 

cipal  from  the  existing  criterion  k  =  Idem,  but  it  reflects  the  latter «s 
physical  significance  more  closely.  It,  like  the  criterion  loddem, 
demands  equality  in  the  values  of  the  index  of  the  adiabatic  curve 
k  in  nature  and  in  the  model. 

In  solving  the  simulation  problem  when  k  /  idem  it  is  natural  to 
turn  to  the  method  of  approximate  simulation.  This  method  recommends 
excluding  from  the  criterion  relation  a  criterion  of  little  effect 
and,  by  means  of  a  practical  check,  confirming  the  suitability  of  the 

Ir 

thus  abbreviated  criteria.  However,  the  criterion  in  the  system 

of  defining  criteria  (1)  is  not  of  little  influence.  Simple  exclusion 
of  it  from  the  standard  relations 


A7™i  tfa 

Gi .  I  '  P»  ’ 


a*  /{  V 
kRT,  ’  k  —  \ )’ 

ua  k  N 
kRT,’  k  —  \J  etc. 


(2) 


leads  to  substantial  errors. 

It  would  seem  Impossible  to  use  the  method  of  approximate  simula¬ 
tion.  The  situation  changes,  however,  if  we  prepare  the  criterional 
relation  In  advance,  which  we  shall  call  generalization  of  the  criteria. 
This  preparation  can  be  done  based  on  the  production  of  criteria  or 
the  quotient  of  their  division  also  being  a  criterion.  Then  the 
criteria  which  it  Is  proposed  to  use  for  approximate  simulation  prove 


capable  of  union  with  the  excluded  criterion,  in  our  case  with  the 

tr 

criterion  The  union  must  be  carried  out  taking  into  considera¬ 

tion  the  influence  of  the  individual  criteria  and  vr.lue  complexes  on 
the  relation  to  be  elaborated,  and  also  observing  the  rules  for  their 
transformation. 

The  most  successful  result  of  using  the  method  of  generalizing 
the  criteria  to  the  relations  (2)  is  obtained  with  the  criterion 
A  .  75%  k 

— — .  In  this  case  the  value  of  i0  =  ART0  is  proportional  to 

Ir 

the  criterion  rr-r  .  Using  the  method  of  generalizing  the  criteria  to 


the  criterion  A  *  ^^Ni,  we  obtain 
'6'io ' 

A  ■  75;V 


A-75N, 


ti  —  I  G  ART q 


A-75N. 

GRT, 


G j  ART „ 

the  value  of  which  does  not  depend  on  the  value  of  k. 

In  conformity  with  the  principle  of  the  theory  of  similarity  and 

dimensionality  the  number  of  criteria  after  transformation  must  not 

A  •  75N, 

change.  Therefore  together  with  the  generalized  criterion _ i, 

k  GiQ 

also  enters  the  standard  equation.  Consequently  the  standard 

equation  will  have  the  form: 


■77>Ni  / pt  ««  k  \ 

GRT,  kRf »’  *~1/' 


(5) 


The  adaptability  of  the 


relationship 

A'75Ni  ...« (Pt 
GRT,  ~~l\p kRT,J 


w 


for  simulating  turbine  work-regimes  on  different  working  substances 
was  confirmed  in  particular  by  treating  the  experimental  data  in  the 
criteria  of  (4)  on  tests  of  an  axial  single-stage  turbine  in  air 
(k  =  1.12),  C02  (k  =  1.27),  and  CCI4  (k  =  1.12)  cited  in  a  certain 
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work  [3].  We  should  note  that  C02  and  CCI4  in  the  tests  had  a  certain 
admixture  of  air.  The  treatment  showed  that  the  turbine  character¬ 
istics  in  the  criteria  of  similarity  in  (4) ,  referring  to  the  turbine 
work-regimes  on  various  working  substances  with  differing  values  of 
the  index  of  the  adiabatic  curve,  coalesce  practically  completely 
when  the  two  determining  criteria  u2/kRTo  and  p2/p0  are  equal  (Pig.  1) • 
Contrariwise,  treating  the  experimental  data  of  the  tests  of  the  same 
turbine  in  air  and  in  C02  in  the  criteria 


etc.  gives  a  10-15$  divergence  in  characteristics  (Pigs.  2,  3).  The 
graphs  confirm  the  adaptability  of  the  standard  relationship  in  (4) 
for  simulation  when  k^idem. 

Thus  the  criterion  relation  in  (4)  may  be  recommended  for  use  in 
simulating  work  regimes  of  turbine  engines  on  different  working  sub¬ 
stances  with  differing  values  of  the  index  of  the  adiabatic  curve. 
Representing  the  turbine  characteristics  in  the  criteria  of  similarity 
in  (4)  enables  us,  from  the  results  of  testing  the  turbine  for  one  of 
the  assimilated  compositions  of  the  working  substance,  to  determine 
the  indices  of  its  work  in  various  set-ups  with  arbitrarily  assigned 
composition  for  the  working  substance  (e.g.,  when  investigating  the 
effect  of  water  injection  in  the  gas  tract  of  gas-turbine  units  on 
the  power  of  the  turbine,  when  developing  atomic  gas-turbine  trans¬ 
portation  units,  etc.),  which  essentially  simplifies  a  number  of 
Investigations  of  power  units  with  turbine  engines. 
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Pig.  1.  Turbine  character¬ 
istics  in  similarity  criteria 

Pj  u*  \ 

Y.'  kRTtj' 

constructed  from  the  results 
of  testing  a  turbine  on  var¬ 
ious  working  substances. 

a  ~  kRT\  ~  °’7;  6  ~kRj \  =  1  >05' 

1)  air;  2)  C02;  3)  CCl* 


A-lSNi  / 
GRT„  -ly 


Pig.  2.  Turbine  character¬ 
istics  in  similarity  criteria 


*-/(*• 


constructed  from  results  of 
testing  a  turbine  on  various 
working  substances  (a,  b  —  see 
Pig.  1).  % 

1)  air;  2)  C02 
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a)  .  b) 

Pig.  3.  Turbine  character¬ 
istics  in  similarity  criteria 

f{  A  .JfM 

\r o'  */?r, ^ * 


/ifir. 


constructed  from  results  from 
testing  a  turbine  on  various 
working  substances  (a,  b  —  see 
Fig.  1).  , 

1)  air J  2)  C02 
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METHODS  AND  CERTAIN  RESULTS  OF  THE  MEASUREMENT 


OF  THE  CRITICAL  LOAD  DURING  THE  TRANSITION 
FROM  FILM  TO  BUBBLE  BOILING 

prof.  B.  S.  Petukhov,  Doctor  of  Technical  Sciences, 
and  S.  A.  Kovalev,  Engineer 
Moscow  Institute  of  Power  Engineering 

The  measurement  of  the  second  critical  loads 
during  the  boiling  of  a  liquid  is  considered,  the 
results  of  an  experimental  study  are  presented, 
and  a  method  of  calculating  the  stability  of  the 
film  and  bubble  regimes  of  boiling  is  proposed. 

In  present-day  steam  generators  and  atomic  reactors,  which  oper 
ate  at  high  thermal  loads,  operating  regimes  accompanied  by  film 
boiling  are  possible.  An  exact  knowledge  is  needed  of  the  so-called 
critical  heat  loads,  at  which  a  change  occurs  in  the  boiling  regime 
of  a  liquid  on  the  heating  surface,  in  order  to  achieve  stable  opera 
tion  of  the  heating  surface  within  given  temperature  limits. 

The  many  investigations  conducted  both  in  the  Soviet  Union  and 
abroad  have  made  it  possible  to  accumulate  a  great  deal  of  experi¬ 
mental  material  concerning  critical  loads  during  the  transition  from 
bubble  to  film  boiling,  qcri . 
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Critical  loads  In  the  transition  from  film  to  bubble  boiling 
have  been  much  less  studied.  Certain  investigations  in  recent  years, 
carried  out  with  the  aid  of  high-speed  motion  pictures,  and  theoreti¬ 
cal  treatments  indicate  that  the  termination  of  film  boiling  occurs 

upon  reaching  the  second  critical  load  q  ,  as  a  result  of  a  loss 

^cr£ 

of  the  stability  of  motion  of  the  vapor  film,  a  stability  related  to 
the  stability  of  the  capillary  waves  at  the  interface  of  the  liquid 
and  vapor  phases.  Certain  experimental  data  on  the  second  critical 
loads  are  also  available.*  However,  these  data  are  very  scanty  and 
show  essential  discrepancies.  Such  a  situation  is,  to  a  certain  extent, 
related  to  the  absence  of  a  sufficiently  developed  method  which  would 
permit  reliable  determination  of  the  second  critical  loads.  On  the 
basis  of  an  analysis  of  the  possible  methods  of  measuring  Qcr2>  this 
article  recommends  one  which  was  developed  at  MEI  (Moscow  Institute 
of  Power  Engineering)  and  which  has  a  number  of  important  advantages. 

The  simplest  and  most  convenient  method  under  laboratory  condi¬ 
tions  is  that  of  electrical  heating,  in  which  a  current  is  passed 
directly  through  the  working  element  of  the  apparatus.  In  the  case 
of  electrical  heating,  film  boiling  may  be  obtained  by  two  methods. 

By  gradually  increasing  the  thermal  load  and  passing  through  the 
first  crisis,  a  transition  to  film  boiling  is  achieved.  After 
attaining  steady  film  boiling  the  load  is  reduced,  until  reverse 
transition  to  the  bubble  regime  of  boiling  takes  place.  This  is  the 
load  taken  as  qcr2.  The  basic  shortcoming  of  this  method  lies  in 
the  fact  that  during  the  transition  to  film  boiling  under  large 

*  M.  V.  Borlshanskiy .  Article  in  the  collection  "Problems  of 
Heat  Transfer  During  a  Change  in  the  Aggregate  State  of  a  Substance." 
Gosenergoizdat,  1955. 
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thermal  loads,  the  temperature  of  the  surface  Increases  sharply. 
Therefore  such  materials  as  nlchrome  and  various  alloy  steels  melt, 
and  their  use  may  thus  lead  to  the  destruction  of  the  surface.  This 
method  requires  the  use  of  refractory  metals  or  graphite. 

However,  through  electrical  heating  it  is  possible  to  obtain  a 

film  regime  of  boiling,  skipping  the  first  crisis,  at  relatively  low 

temperatures  and  heat  loads  on  the  heating  surface.  To  do  this,  the 

heating  surface  is  warmed  up  by  an  electric  current  in  the  vapor  phase 

to  a  temperature  higher  than  the  temperature  of  the  second  boiling 

crisis  t  =  t_  +  t  and  is  then  lowered  into  the  liquid  phase  under 
s  cr  <=- 

a  thermal  load  somewhat  greater  than  that  of  the  second  critical 
load.  For  example,  in  order  to  obtain  film  boiling  in  water  at 
atmospheric  pressure,  it  is  sufficient  to  heat  the  surface  to  only 
550  -  400°  C. 

The  method  of  electrical  heating  with  transition  through  the 
first  critical  load  was  used  in  the  work  of  V.  M.  Borishanskiy  in  an 
investigation  of  the  second  critical  loads  of  water,  isooctane,  and 
benzene  under  conditions  of  free  convection  at  pressures  from  i  to 
16  atm  abs.  Graphite  rods  with  diameters  of  2.2  mm,  secured  in  copper 
clamps,  were  used  as  the  experimental  heating  surfaces.  In  these 
experiments  it  was  established,  in  particular,  that  the  second  criti¬ 
cal  load  for  water  at  atmospheric  pressure  is  2  •  105  kcal/m2  •  hr. 

It  should  be  noted  that  in  the  experiments  of  V.  M.  Borishanskiy 
there  Is  no  doubt  that  significant  leakages  of  heat  through  the  ends 
of  the  graphite  rods  took  place,  and  these  might  have  substantially 
distorted  the  values  of  the  second  critical  load.  As  an  example, 
let  us  consider  a  sufficiently  long  rod,  the  ends  of  which  are  secured 
In  massive  copper  current  carriers,  while  on  the  surface,  film  boiling 
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is  maintained  by  electrical  heating.  Since  the  ends  of  the  rod  will 
have  a  ]ow  temperature  corresponding  to  the  temperature  of  the  massive 
clamps,  at  these  places  on  the  surface  of  the  rod  there  must  be  a 
bubble  regime  of  boiling.  Thus,  near  the  ends  of  the  rod  there  will 
be  zones  with  two  regimes  of  boiling  —  bubble  and  film  —  with  very 
different  heat-transfer  coefficients  and  sharply  differentiated 
surface  temperatures.  As  a  result  of  this  simultaneous  occurrence 
of  two  regimes  of  boiling  at  the  ends  of  the  rod,  with  a  reduction 
of  the  load  ,  the  film  starts  to  recede  to  the  center  of  the  rod 
long  before  the  attainment  of  the  second  critical  density  of  the  heat 
flux  q  ,  on  account  of  a  reduction  in  the  temperature  of  the  heating 

OX  2 

surface.  Regardless  of  the  length  of  the  rod  and  the  relation 
between  the  quantity  of  heat  given  off  from  the  surface  and  the  end 
losses,  this  process  will  continue  until  the  film  is  removed  from  the 
entire  surface  of  the  rod.  Under  these  conditions  (on  account  of  a 
reduction  in  the  temperature  by  axial  heat  flows)  the  local  heat  loads 
from  the  surface  in  the  zone  of  film  removal  will  be  significantly 
less  than  the  heat  loads  obtained  in  the  experiments  and  calculated 
according  to  the  electrical  power  input  over  the  entire  surface  of  the 
rod.  Consequently,  on  a  sufficiently  long  rod  it  is  possible  to 
study  heat  transfer  even  in  the  presence  of  end  losses,  but  it  is  not 
possible  to  study  critical  thermal  loads. 

For  these  reasons  the  data  obtained  by  V.  M.  Borishanskiy  and 
other  authors  by  the  methods  described  can  hardly  be  considered 
trustworthy.  They  undoubtedly  require  verification.  For  this  purpose 
special  experiments  were  carried  out  by  MEI. 

The  experiments  were  carried  out  with  distilled  water  at  atmos¬ 
pheric  pressure  under  conditions  of  free  convection.  Horizontal  wires 
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and  tubes,  through  which  an  electric  current  was  passed,  served 
as  the  heating  surface.  In  order  to  obtain  film  boiling,  the  above- 
described  method  of  heating  the  specimen  to  the  vapor  phase  before¬ 
hand  was  used. 


Special  consideration  was  given  to  obtaining  a  constant  tempera 
ture  field  throughout  the  length  of  the  specimen.  Since  it  was  not 
possible  to  secure  the  wires  or  tubes  in  such  a  way  that  there 


would  be  no  heat  leakages  in  the  current  conductors,  the  copper 
clamps  were  withdrawn  into  the  vapor  phase  (Pig.  1).  This  made  it 
possible  to  eliminate  the  temperature  irregularity  along  the  part  of 


the  heating  element  located  in  the  liquid  phase  by  selecting  the 
height  of  the  protruding  ends . 


Pig.  1.  The  temperature  dis¬ 
tribution  over  the  length  of 
a  stainless-steel  specimen  with 
a  diameter  of  2  mm  at  q  = 

=  4o  •  103  kcal/m2  •  hr. 

1  —  1,  —  !  In  .tin;  2— 1,-1  -  .tin;  3-1,— l— \  mm. 


In  order  to  calculate  the 
temperature  distribution  along  the 
length  of  the  specimen,  we  shall 
assume  that  the  temperature  does 
not  vary  over  the  cross  section 
of  the  specimen.  Then  the  differ¬ 
ential  equation  of  the  temperature 
field,  which  is  easily  obtained 
from  the  heat  balance,  has  the 
form: 


_ 0 ,  SO  l -ft  nu 

'  H  A  ’ 

_o,sn/^  »«/, 

1  /2a  'i  71  ’ 


(1) 


where  t  =  temperature  of  the  wire  or  tube  at  a  given  point  (°c); 

t  =  saturation  temperature  (°C); 
s 

p  =  a  +  bt  =  resistivity  (ohm.m); 

I  =  current  (amps); 


a  =  heat-transfer  coefficient  (kcal/m2  •  hr  •  °C); 

X  =  coefficient  of  thermal  conductivity  (kcal/m  •  hr  •  °C); 
f  =  cross-sectional  area  of  the  specimen  (m2); 
u  =  perimeter  of  the  specimen  (m)j 

Since  the  constants  3  and  v  contain  the  heat-transfer  coefficient, 
they  will  have  different  numerical  values  for  elements  of  the  surface 
located  within  the  liquid  and  vapor  phases.  Indeed,  for  elements 
of  the  surface  located  in  the  liquid,  where  the  heat-transfer  coeffi¬ 
cient  is  high  even  for  film  boiling,  3  <  0,  while  for  elements  of  the 
surface  located  in  the  vapor  phase,  where  a  is  small,  usually  3  >  0. 
Accordingly,  the  solutions  of  Eq  (1)  will  have  the  following  form: 

for  elements  of  the  surface  of  the  specimen  located  in  the 
liquid  (3  <  0)  i.e.,  over  the  interval  0  <  x  <  i. 


t 
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‘shi'l/  -ru*chV~u  +  P  ’ 
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sli  V  tu  r  pi 

*  eh  ■  oosh  and  sh  ■  sinh 

for  elements  of  the  surface  of  the  specimen  located  In  the  vapor 


phase  (3  >  0)  i.e.,  over  the  interval  i  <  x  <  u 


t = c  sin  y'? i(-v'— o  i  q  cos  v p,  (a- — i) 
1  sin —  ‘  ‘ 


Y  i 


cos /[,(/,  —  !)  ?i 
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Here  the  total  length  of  the  specimen  is  taken  equal  to  2llt  while 
the  length  of  the  section  of  the  specimen  immersed  in  the  liquid  is 
2l . 


In  the  particular  case  where  3=0  the  variation  in  the  tempera¬ 
ture  will  be  parabolic. 

In  order  to  find  the  constants  of  integration,  it  is  necessary 
to  assign  boundary  conditions. 

Since  the  temperature  distribution  must  be  symmetrical,  at  the 
center  of  the  specimen  (i.e.,  for  x  =  0)  dt/dx  =  0  (Fig.  1). 

The  temperature  of  the  ends  of  the  specimen  is  equal  to  the 
temperature  of  the  clamps  tx.  In  other  words,  for  x  =  ix,  t  =■  tx. 


It  may  be  easily  measured  experimentally  and  is  usually  close  to  the 
saturation  temperature. 

In  order  to  obtain  a  description  of  the  temperature  field  over 
the  entire  length  of  the  specimen,  it  is  necessary  to  reconcile  the 
solutions  to  Eqs.  (2)  and  (3),  i.e.,  for  x  =  i  the  following  conditions 
must  be  satisfied: 


K 


-  /  —  I  ( x>t  I  x=r 


dtr 
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dx 
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lx>! 


dx 
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After  simple  transformations  we  arrive  at  the  values  of  the 
constants  in  Eqs.  (2)  and  (3): 


D,  =  0; 


V,  ,  Y,  ^  1 

1 

V  1  1  /  cos  /?,#,-/) 

Lvsl 

!«*• 

C,  =  [(/,+ k)  th  /p7  •  tg  v%  (/,  -  0  -  ^ 

A; 

c3=[/,+|f+(j+^)/fth/F/X 

XsinkFTk-0]  A,  . 


(2a) 


(5a) 


^  =  [l  +  j/£tli/p7tg  /&(/,  -/)]“. 

(note:  th  ■  tanh  and  tg  *tan) 

Analysis  of  Eqs.  (2)  and  (2a)  indicates  that  it  is  possible  to 
select  a  length  of  the  protruding  ends  ( h  -  l)Q,  such  that  the 
second  term  on  the  right-hand  side  of  Eq,  (2)  goes  to  0  and  the 
temperature  on  the  entire  section  of  the  specimen  immersed  in  the 
liquid  will  be  the  same.  Assuming  Ds  =  0,  from  Eq.  (2a),  we  find: 


It  is  clear  from  simple  physical  considerations  that  if  the 
section  of  the  specimen  protruding  from  the  liquid  has  a  length 
greater  than  (li  -  l )Q,  then  the  specimen  in  the  vapor  phase  will  be 
superheated,  and  if  it  is  smaller  than  -  l  )Q,  then  the  section 
of  the  specimen  located  in  the  liquid  will  be  cooled  as  a  result  of 
heat  flowing  off  into  the  clamps.  The  same  results  follow  from  Eqs. 
(3)  and  (3a).  As  an  example.  Pig.  1  shows  the  temperature  distribu¬ 
tion  along  a  two-millimeter  stainless-steel  wire,  calculated  accord¬ 
ing  to  Eqs.  (2a)  and  (3a)  for  three  values  of  u  -  i . 

The  graphs  shown  in  Pig.  1  clearly  indicate  that  the  relatively 
small  deviations  of  h  -  i  from  (ix  -  i )  significantly  change  the 
temperature  field  of  the  specimen  in  the  vapor  phase  and  on  a  small 
section  of  the  specimen  in  the  liquid  near  the  free  level,  but  have 
practically  no  affect  on  the  entire  central  part  of  the  specimen. 
Therefore,  experiments  on  the  investigation  of  q  may  be  carried 
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out  with  protruding  ends  longer  than  (lx  -  l)Q.  In  this  case 
removal  of  the  film  will  always  begin  in  the  central  part  of  the 
specimen. 

A  profile  with  a  monotonic  reduction  of  temperature  toward  the 
ends  for  lx  -  i  <  (ix  -  i )Q  as  was  demonstrated  above,  is  completely 
unsuitable  for  a  study  of  q  ,  but  a  constant  profile  for  lx  -  l  = 

vX  2 

=  (z x  -  l )Q  is  difficult  to  realize.  If  in  the  process  of  the 
experiments  it  should  turn  out  that  the  calculated  height  (u  -  j ) 
was  underestimated  and  the  removal  of  the  film  begins  at  the  points 
where  the  specimen  emerges  from  the  liquid,  then  it  may  be  easily 
detected  visually  and  corrected  by  changing  the  height  of  the  protrud¬ 
ing  ends. 
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Figure  2  shows  a  schematic  of  the  apparatus.  It  consists  of 
a  stainless-steel  tank  which  has  a  window  in  the  front  wall  for  visual 
observations.  The  experimental  specimen  is  secured  in  copper  clamps 
which  are  mounted  on  a  textolite  panel  and  insulated  from  the  steel 
cover  of  the  tank  by  textolite  bushings.  During  the  experiment  the 
panel  can  easily  be  moved  upwards;  then  the  specimen  is  in  the  vapor 
phase.  The  water  in  the  apparatus  is  maintained  at  the  saturation 
temperature  by  an  electric  heater  located  beneath  the  bottom  of  the 
tank.  In  order  to  avoid  the  Influence  of  convection  currents  of  the 
liquid  on  the  vapor  film  surrounding  the  specimen  there  is  a  special 
shield  between  it  and  the  bottom.  In  order  to  measure  the  voltage 
drop  over  the  experimental  section  of  the  specimen,  two  potentiometer 
probes  0.2  mm  in  diameter  were  welded  to  its  upper  generatrix  at  a 
distance  of  about  100  mm  from  each  other.  Specimens  with  diameters 
of  0.5  and  1.4  mm  were  prepared  from  constantan  wire,  so  as  not  to 
distort  the  temperature  field  with  the  potentiometer  wires.  In  this 
case  the  thermal  loads  were  determined  from  the  current  and  the 
resistance. 

The  experimental  procedure  was  as  follows.  After  the  water  was 
brought  to  a  boll,  the  specimen  was  raised  to  the  vapor  phase  and  was 
heated  by  an  electric  current  to  a  temperature  of  550-400°.  The 
heated  specimen  with  a  thermal  load  greater  than  the  second  critical 
load  was  lowered  into  the  water,  and  film  boiling  set  in  on  its  upper 
surface.  The  thermal  load  was  reduced  by  small  stages,  the  magnitude 
of  which  did  not  exceed  1%  of  the  load.  After  each  reduction  of  the 
load  and  attainment  of  a  steady-state  regime,  measurements  were  made 
of  the  current  and  the  voltage  drop  over  the  experimental  section. 

The  regime  at  which  removal  of  the  film  sets  in  was  considered  to 


be  the  critical  regime  corresponding  to  the  second  critical  load. 

Such  regimes  were  often  repeated.  Control  experiments  were  conducted 
systematically,  in  which  film-boiling  regimes  immediately  preceding 
the  crisis  were  maintained  for  10-20  minutes.  These  experiments 
served  as  good  confirmation  of  the  reliability  of  the  determination 


Pig.  2.  Experimental  apparatus;  i)  tank, 
s  =  380  x  190  mm;  2)  experimental  speci¬ 
men;  5)  potentiometer  leads;  4)  current 
carriers;  5)  Insulating  bushings;  6) 
tank  cover;  j)  textolite  panel;  8) 
electric  heater;  9)  protective  shield; 

10)  window. 


In  each  experiment  the  location  of  film  removal  was  monitored 
visually.  Usually  the  film  began  to  come  off  at  an  arbitrary  point 
of  the  horizontal  part.  The  beginning  of  film  removal  from  the  ends 
of  the  specimen  attested  to  the  presence  of  heat  leakages  along  the 
specimen,  and  such  experiments  were  rejected. 

During  the  experiments  the  following  picture  of  film  removal 
was  observed.  Prom  thin  wires  the  film  comes  off  rapidly  (fractions 
of  seconds)  with  a  characteristic  crackling.  The  film  comes  off 


slowly  from  thick  stainless-steel  pipes  with  a  great  heat  capacity 
and  we  may  readily  observe  how  at  critical  loads  the  film-boiling 
regime  is  replaced  by  unstable  film  oscillations,  which  correspond 
to  a  transition  region  and  which  then  give  way  to  the  bubble-boiling 
regime . 

The  experimental  results  are  presented  in  Fig.  3  in  the  form  of 
a  dependence  of  the  second  critical  thermal  load  on  the  diameter  of 
the  specimen  (in  logarithmic  coordinates).  A  study  of  the  experimental 
data  shows  that  the  second  critical  load  decreases  with  an  increase 
in  the  diameter  of  the  specimen.  A  dependence  of  the  critical  loads 
on  the  material  of  the  surface  was  not  revealed.  The  measurement 
results  are  presented  in  the  table. 
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Fig.  3>  Dependence  of  the  critical  load 

on  the  diameter  of  the  specimen:  l) 

q  „  ,  data  of  the  authors:  2)  data  of 
^cr2 

V.  M.  Borishanskiy;  3)  qequl  experi¬ 
mental;  4)  q.equl  calculated. 

The  values  of  q  obtained  make  it  possible  to  determine  by 
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calculation  the  thermal-load  values  of  practical  interest,  at  which 
the  bubble  and  film  regimes  of  boiling  will  be  in  equilibrium  on  a 


wire  or  tube.  Let  us  call  this  load  the  equilibrium  load  (^equi^* 

For  heat  flows  greater  than  QeqUj_  film-boiling  regimes  which  arise 
accidentally  will  propagate  over  the  entire  surface  of  the  specimen. 
For  loads  less  than  QeqUj_  &nd  with  the  simultaneous  existence  of  film 
and  bubble  boiling  on  the  surface  of  the  specimen,  film  removal  occurs 
as  a  result  of  axial  thermal  outflows. 

In  order  to  calculate  the  equilibrium  lead,  let  us  use  the  follow 
ing  model.  Developed  bubble  boiling  is  maintained  on  the  outer  surfac 
of  a  tube  of  length  21.  Let  us  find  the  thermal  loads  at  which  both 
regimes  of  boiling  will  be  stable. 
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Fig.  4.  Temperature  distribution  over 
the  length  of  a  specimen  under  an 
equilibrium  load. 


It  is  known  that  in  the  case  of  film  boiling  the  dependence  of 
the  heat  transfer  on  the  load  is  slight.  Therefore  it  may  be  assumed 
that  in  the  region  of  film  boiling  the  heat-transfer  coefficients  are 
constant  over  the  length  of  the  tube  and  are  equal  to  As  for 

the  transition  region  and  the  region  of  bubble  boiling,  the  heat- 
transfer  coefficients  may  vary  over  the  length  of  the  tube.  However, 
for  simplicity  of  calculation,  we  shall  take  the  heat-transfer  coef¬ 
ficients  in  both  of  these  regions  as  being  equal  to  the  heat-transfer 
coefficients  for  ordinary  bubble  boiling  abUbble  at  a  given  load. 

Since  the  heat-transfer  coefficients  for  the  film  and  bubble  regimes 
of  boiling  are  different,  the  temperature  field  of  the  specimen  is 
described  by  two  equations  of  type  (2),  which  must  be  reconciled  at 
the  point  of  change  in  the  boiling  regime  with  respect  to  temperatures 
and  temperature  gradients.  As  the  boundary  conditions  we  shall 
assume  that  the  temperature  gradients  in  the  center  of  the  tube  and 
in  the  zone  of  bubble  boiling  far  from  the  interface  are  equal  to 
zero.  We  shall  also  consider  as  given  the  film -removal  temperature, 
which  is  determined  by  the  experimental  values  of  the  second  critical 

load  t  =  t  +  . 

c  si ilm 

The  temperature  distribution  in  the  region  of  film  removal,  thus 
calculated,  is  shown  in  Pig.  4.  Prom  the  general  solution  we  obtain 


the  following  equation; 


p+4-  V  thVjfi, 

^  _  n  Y  y  Pi _ _  * 


*cth  *  coth  and  th  ■  tanh 


where 


’  the  temperature  of  the  specimen  in  the  region  of 
?  film  boiling  in  the  absence  of  axial  heat  flows; 

ir~ W>bie  *  the  temperature  of  the  specimen  in  the  region  of 

e  bubble  boiling  in  the  absence  of  axial  heat  flows. 
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For  a  sufficiently  long  tube  or  wire,  the  equation  becomes  simply: 

k  —  hum  1 


'bubEli  filml  +  ]/"  JL 

f  ft 


(6) 


If  the  electrical  resistance  of  the  tube  does  not  depend  on  the  temper¬ 
ature,  Eq.  (6)  assumes  the  form: 

i  ~  _ ! _ 

“ 


iubble^  film  1  + 


‘film 

‘bubble 


(7) 


Using  Eq.  (j)  and  the  relationships 


qequi  ~  “film— film 


-  t)  =  ay 


(t> 


-  t  )  and  a. 


s '  "hubble-bubble  “s'  “““  “"bubble  " 

=  f(tbubble  -  tg)  it  is  easy  to  obtain  tfllm  and  the  equilibrium  load. 


For  a  two-millimeter  wire  with  a,  , ,  , „  =  2 

bubble 


10*  and  cfllm  =  'too 


kcal/m2  •  hr  •  °C  (afjLim  taken  according  to  data  of  a  specially 
conducted  experiment)  we  obtain  from  Eq.  (7)  Qequi  -  25^ "103  kcal/m2* hr. 

The  results  of  the  calculation  were  verified  experimentally.  The 
experiments  were  conducted  on  a  horizontally  stretched  nlchrome  wire 
250  mm  long  and  2  mm  in  diameter.  In  order  to  obtain  film  boiling, 
the  wire  was  heated  in  the  vapor  phase  and  then  plunged  into  the 
liquid  phase.  The  load  was  gradually  reduced,  until  the  film  began 
to  come  off  the  ends  of  the  wire.  When  the  film  covered  about  l/j5 
of  the  length  of  the  wire,  a  load  was  selected  at  which  film  removal 
ceased  and  the  interface  between  the  regimes  of  boiling  became  station¬ 
ary.  This  load  was  taken  as  Qequ;j_*  When  the  load  was  increased,  the 
film  slowly  spread  out  over  the  entire  surface,  while  with  a  decrease 
in  the  load  further  removal  of  the  film  occurred.  The  experimental 
results  (<lequi  =  276 *103  kcal/m2 *hr)  corroborate  the  calculated  value 
of  the  equilibrium  load. 

Thus,  in  order  to  calculate  equilibrium  thermal  loads,  it  is 
sufficient  to  know  the  exact  values  of  the  heat-transfer  coefficients 
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and  the  second  critical  loads.  Analogous  calculations  of  equilibrium 
loads  may  be  made  for  a  heating  surface  of  any  other  shape. 

The  experimental  data  obtained  by  V.  M.  Borishanskiy,  as  follows 
from  Pig.  5  and  also  from  the  nature  of  the  method  employed,  correr 
sponds  more  to  ai  equilibrium  thermal  load  than  to  the  second  boiling 
crisis . 


Conclusions 

1.  A  method  has  been  developed  for  the  measurement  of  the  second 
critical  loads  during  the  boiling  of  a  liquid  under  conditions  of 
free  convection,  a  method  which  makes  it  possible  to  obtain  reliable 
experimental  data. 

2.  Measurements  were  made  of  the  second  critical  load  on  hori¬ 
zontal  tubes  and  wires  with  diameters  ranging  from  0.5  to  10  mm  for 
water  under  conditions  of  free  convection  in  a  large  volume.  It  is 
shown  that  the  published  data  concerning  second  critical  loads  exceed 
the  actual  values  several  times  over. 

5.  It  is  shown  that  the  second  critical  load  deereases  with  an 
Increase  In  the  diameter  of  the  wire  or  tube. 

4.  The  material  of  the  heating  surface  was  not  observed  to  have 
any  effect  on  qcr2. 

5.  It  was  shown  that  the  experimental  values  of  qcr>2  obtained 
make  it  possible  to  calculate  the  equilibrium  loads  which  constitute 
the  upper  limit  of  stable  bubble  boiling  and  the  lower  limit  of  stable 
film  boiling  on  a  given  surface. 
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METHODS  OP  GENERALIZING  EXPERIMENTAL  DATA 


ON  CONVECTIVE  HEAT  TRANSFER  DURING  THE 
MOVEMENT  OP  GAS  IN  THE  INITIAL 
REGION  OP  A  CONDUIT 


Prof.  S.  I.  Kosterin,  A.  I.  Leont’ yev, 
and  V.  K.  Fedorov 


A  new  method  for  generalizing  the  experimental 
data  on  convective  heat  transfer  in  the  turbulent 
boundary  layer  of  a  gas  is  proposed.  A  universal 
law  of  heat  transfer  for  the  turbulent  boundary 
layer  of  a  compressed  gas  is  derived. 


Research  on  the  processes  of  convective  heat  transfer  and  of 
resistance  in  the  motion  of  a  gas  in  relatively  short  cylindrical 
conduits  is  at  present  timely  enough  from  both  the  practical  and  the 
scientific  point  of  view.  There  are  sufficiently  reliable  analytical 
solutions  for  the  case  of  laminar  boundary  layer  in  the  initial 
region  of  a  cylindrical  conduit.  Methods  of  calculating  this  layer, 
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however,  in  elaborating  which  experimental  investigations  exert  great 
influence,  are  of  great  practical  interest.  We  cannot  consider  as 
successful  the  methods  for  this  sort  of  generalization  recommended 
in  the  literature  [1-3]. 

One  article  [l]  cites  interesting  experimental  data  on  heat 
transfer  and  on  resistance  in  cases  of  sufficiently  Intense  heating 
of  the  gas.  The  method  of  generalization  adopted  by  the  authors, 
however,  is  of  a  very  provisory  nature;  and  the  proposed  criterional 
equations  have  a  limited  area  of  application.  Proceeding  from 
elementary  considerations  we  may  show  that  the  criterional  equation 
of  heat  transfer  in  the  motion  of  a  compressed  gas  in  the  initial  sec¬ 
tion  of  a  cylindrical  conduit  has  the  form 

Nu  =  f  (Re,  Pr,  x,  M,  \)  .  ( l) 

The  authors  [1]  write  Eq.  1  in  the  form 

Nu  =  c  RemPr°*4(  )  p,  (2) 

V 

where  the  coefficients  c,  m,  and  p  for  the  initial  section  are  func¬ 
tions  of  x/D  and  these  functions  are  derived  by  experimentation  and 
presented  in  the  form  of  separate  graphs.  It  is  quite  clear  that  the 
proposed  method  of  calculation  will  give  correct  results  only  for  those 
heat-supply  relationships  which  occurred  in  the  experiments  of  the 
authors  of  the  article  [1].  We  should  also  note  that  under  the  common 
influence  of  several  criteria  it  is  difficult  to  fix  the  effect  of  each 
criterion  separately;  and  the  nature  of  spread  of  the  experimental 
points  near  the  path  of  the  criterional  equation  may  not  be  used  for 
final  evaluation  of  the  criterional  equation's  correctness,  since 
under  these  conditions  mutual  compensation  of  errors  is  possible.  The 


one-dimensional  flow  chart  for  the  gas  in  the  initial  section  of  the 
conduit  adopted  in  the  article  [1]  does  not  correspond  to  the  actual 
picture  of  the  process  and  can  lead  to  incorrect  conclusions.  In 
particular,  the  very  hazardous  conclusion  that,  when  isothermal  con¬ 
ditions  are  lacking,  the  basic  relationship  of  the  hydrodynamic  theory 
of  heat  transfer  fails  is  the  result  of  treating  the  experimental 
data  on  the  coefficient  of  friction  in  accordance  with  this  one¬ 
dimensional  flow  chart.  If  we  treat  the  data  in  accordance  with  the 
two-dimensional  flow  chart,  no  failure  of  the  similarity  is  observed. 

We  should  also  note  that  the  proposed  crlterional  equation  is  true  only 
for  those  conditions  at  the  conduit  entry  which  occurred  in  the  experi¬ 
ments  in  the  article  [1]. 

An  interesting  method  for  generalizing  the  experimental  data  on 
heat  transfer  in  the  initial  section  of  a  cylindrical  conduit  is  pro¬ 
posed  in  another  article  [2].  This  method  is  convenient  for  treating 
the  experimental  data  and  the  practical  calculations  when  a  law  of  heat 
supply  along  the  length  of  the  conduit  is  given.  The  principal  defect 
of  this  method  is  the  necessity  of  deriving  its  crlterional  equation 
for  every  heat-supply  relationship. 

Experimental  research  on  heat  exchange  in  conduits  at  supersonic 
gas  velocities  when  entering  the  conduit  Is  of  great  interest  In  other 
articles  [J>,  4].  In  this  case  we  must  take  into  consideration  the 
effect  of  gas  compressibility  on  the  coefficients  of  heat  transfer  and 
this  substantially  complicates  crlterional  treatment  of  the  experimental 
data.  We  should  note  that  at  supersonic  velocities  of  the  gas  on  en¬ 
trance  into  the  conduit  we  are  unable  to  obtain  supersonic  velocities 
of  gas  flow  for  a  great  extent  and  therefore  it  is  almost  always  nec¬ 
essary  to  to  deal  with  the  region  of  hydrodynamic  and  thermal  stabiliza¬ 
tion.  Besides  that,  there  arise  essential  difficulties  connected  with 
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The  dependence  of  St0  on  Peg for  the  initial 
section  of  a  conduit.  Data  from  l)Lel!chuk 
and  Ity-adyakin  [1];  2)  Petukhov,  Detlav  and 

Kirillov  [11]J  5)  Bradfleld  and  Decoursin 

[10],  Law  of  heat  transfer  for  i)  turbulent 
boundary  layer  ^  0.014  ; 

st°  ‘  ^fs^s:  5 


II)  laminar  boundary  layer  Sto  = 


0,22 
Pe0  Pr1^ 


obtaining  an  undisturbed  supersonic  stream  entering  the  conduit;  there¬ 
fore  the  accuracy  of  the  experiments  under  these  conditions  is  sub¬ 
stantially  lower  than  in  subsonic  velocities  of  gas  motion.  In  the 
article  [4]  the  bulk  of  the  experimental  points  are  distributed  in  the 
transitional  region;  therefore  the  authors  have  renounced  criterional 
treatment  of  the  experimental  data. 

One  article  [3]  cites  the  extensive  experimental  material  on  the 
coefficients  of  heat  transfer  at  supersonic  velocities  of  gas  motion 
in  a  conduit.  Here  the  authors  attempt  to  generalize  the  experimental 
data  in  the  form  of  criterional  equations  based  on  one-  and  two-dimen¬ 
sional  simulations  of  flow  in  the  conduit.  Under  these  conditions  the 
number  of  identifying  criteria  mutually  influencing  the  local  Nusselt 
criteria  increases  to  four  (ReD,  Xj  Tw,  M  )  .  Experiments  have  been 
conducted  under  comparatively  slight  intensities  of  heat  transfer; 
therefore  the  criteria  Re^,  x,  and  M  exert  a  basic  influence  on  the 
change  in  the  coefficients  of  heat  transfer  along  the  length  of  the 
conduit;  and  in  this  the  distribution  of  M  along  the  length  is  a  func¬ 
tion  of  ReD*  It  Is  extremely  hazardous  to  say  that  under  these  complex 
conditions  the  final  criterion  equation  takes  into  consideration  the 
effect  of  every  criterion  separately.  In  particular,  the  authors  of 
the  article  [3],  when  determining  the  influence  of  compressibility, 
assume  in  essence  that  the  effect  of  criterion  x  on  the  local  values 
of  the  Nusselt  criteria  is  the  same  at  both  subsonic  and  supersonic 
velocities,  which  Is  not  in  accord  with  the  realities  of  the  situation, 
since  the  development  of  the  boundary  layer  along  the  length  of  the 
conduit  is  essentially  a  function  of  M  on  entry. 

Thus  we  may  infer  that  at  present  we  have  no  reliable  and  physi¬ 
cally  based  method  of  generalizing  experimental  data  on  convective 


heat  transfer  in  gas  motion  in  the  initial  section  of  a  conduit. 

It  is  quite  natural  to  use  the  theory  of  local  simulation  in 

investigating  these  heat-transfer  processes.  Two  works  [ 5 ,  6]  expound 
the  basic  ideas  of  this  theory.  According  to  this  theory  the  aim  of 
the  experiment  is  to  establish  the  laws  of  heat  transfer  and  resistance 
in  the  turbulent  boundary  layer;  but  the  effect  of  the  different 
external  conditions  (distribution  of  pressures,  wall  temperatures, 
thermal  currents,  etc.)  are  accounted  for  by  the  boundary-layer 
equations.  The  equation  for  the  thermal  boundary  layer  for  the  move¬ 
ment  of  a  gas  in  the  initial  section  of  a  conduit  may  be  written  in 
the  following  form: 


(.ivate  * . st  i v  • 
ax  MV 
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P  « 

w*.  V 


W  ~*W  jfc  -fc  ) 


_ l  w _ .  pe 

’  D  a« 
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where  Too  is  the  stagnation  temperature  in  the  center  of  the  stream; 

D)  diameter  of  the  conduit;  T*  )  equilibrium  temperature  of  the  wall; 

T*  )  stagnation  temperature  in  the  boundary  layer;  Tw  )  wall  temperature; 

Po,  u0  )  density  and  velocity  of  the  gas  in  the  center  of  the  stream; 

q  )  thermal  current  at  the  wall;  a0o*  *-oo  )  coefficients  of  temperature 
w 

conductivity  and  heat  conductivity  at  stagnation  temperature  T0q»  p,  u) 
density  and  velocity  of  the  gas  in  boundary  layer. 
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Prom  Eq.  1  we  get 


or 
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Pc„ 


Qx 


(3) 


where  , 

Qx=$9DfydLc- ' 
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We  must  find  the  parameters  of  the  gas  in  the  center  of  the 
stream  in  order  to  determine  by  experiment  the  local  values  of 
Stanton's  criterion.  If  the  static  pressures  along  the  conduit  are 
measured,  the  parameters  of  the  gas  in  the  center  of  the  current  are 
immediately  determined  from  the  gas-dynamic  functions.  But  if  only 
thermal  measurements  are  made  in  the  experiments,  we  may  then,  with 
sufficient  accuracy  for  the  problem  posed,  take,  as  do  Leont'yev  and 
Fedorov  [ 7] , 


Pe  = 
Pe5  1 


:  Pr  Re*; 

_  Po‘‘ 9^ 


P(|» 


(*) 


where  $  is  the  size  of  the  loss  of  momentum.  Then,  using  the  dis¬ 
continuity  equation,  we  obtain  (loc.  cit.) : 


where 


UhH  —  Ren + 4/iReft , 
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and 
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In  subsonic  and  perisonic  gas  flow  h  =  l.JT,  .  which  is  known 
from  the  literature  [5>  6,8].  Thus,  based  on  the  usual  measurements 
made  in  experimental  research  on  heat  transfer  we  may  from  Eqs.  2  and 
6,  determine  the  local  values  of  the  Stanton  and  Pe'clet  criteria  and 
establish  the  law  of  heat  exchange.  The  figure  shows  the  results  ob¬ 
tained  from  treating  the  experimental  data  on  heat  transfer  in  the  ini¬ 
tial  section  of  a  conduit  on  a  plate  and  a  cone  in  accordance  with  the 
proposed  methods.  The  effect  of  anisothermicity  and  compressibility 
on  the  local  values  of  Stanton's  criteria  is  accounted  for  in  accord 
with  the  equation  suggested  by  Kutateladze  and  Leont 'yev [9] .  All  the 
experimental  points  obtained  for  the  various  laws  of  heat  supply  in  the 
different  experimental  set-ups  have  been  placed  on  one  curve,  which 
may  be  described  by  the  following  equation: 


where 


V  *CpT<>0 


St 


0,014(1 


(7) 


Utilizing  Eqs.  5  and  7  and  integrating  Eq.  1  we  may  obtain  the 
solution  for  the  cases  where  Tw  =  const  and  the  law  of  heat  supply 
is  given,  cp.  Leont'yev  and  Fedorov  [7]. 

The  advantage  of  the  proposed  methods  of  generalization  of  the 
experimental  data  versus  criterion  treatment  is  our  ability  to  ex¬ 
clude  the  effect  which  criterion  x,  basically  determining  the  distri¬ 
bution  of  the  coefficients  of  heat  transfer  along  the  length  of  the  con¬ 
duit,  exerts  on  the  laws  of  heat  transfer  and  to  ascertain  the  effect 
of  anisothermicity  and  compressibility  in  their  pure  form.  In  addi¬ 
tion,  the  proposed  law  of  heat  transfer  has  a  universal  nature  and  all 
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the  multiplicity  of  boundary  conditions,  i.e.,  the  direct  effect  of 
x  and  of  the  law  of  heat  supply  on  the  distribution  of  local  coeffi 
cients  of  heat  transfer,  are  taken  care  of  by  the  equation  of  the 
thermal  boundary  layer. 
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